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ABSTRACT 

We show that the first structures that form in the universe should spontaneously 
generate magnetic fields. No primordial seed field is required for this "first dynamo" . 

Although the first dynamo starts with kinetic plasma instabilities, we argue that 
an adequate magnetohydro dynamic description might be possible via a simple trick. 
This should allow a numerical study of the effect of the first magnetic fields on the first 
baryonic objects. 

1. Introduction 

It has already been proposed that the very first objects that form in the universe should 
generate magnetic fields (Gruzinov 2001, Rees 2006). Here we offer a concrete scenario. In our 
scenario the first magnetic fields appear before the first stars, and the dark age ends with the 
synchrotron light from the first shocks. This light is probably unobservable. But the first magnetic 
fields can, in principle, change the first star or black hole formation. 

Consider a naive version of the magnetic field evolution at high redshift, z ~ 100, when the 
universe contains no stars or quasars, and the very first baryonic structures begin to form. The 
gas turbulence generated by the structure formation amplifies the magnetic field, because the gas 
is a good conductor. The growth time of the field is of order the dynamical time, and the first 
structures amplify the primordial magnetic field by a factor of a few. Since the primordial field 
must be very weak0, we must conclude that the first structure formation is virtually unaffected by 
the magnetic field. 

This logic is flawed, because the exponentiation time of the magnetic energy is actually much 
shorter than the dynamical time. We explain the main idea in §2, give the full scenario in §3, where 
we also provide the relevant numerical estimates. In §4 we propose a simple method which should 
allow to use numerical magnetohydrodynamics to study the first magnetic fields. 



^ For all we know, there might be no primordial magnetic field at all. There are some speculative theories of 
primordial magnetogenesis, but they too predict only a weak field. 



- 2 - 



2. The main idea 

Our full scenario is too complex and multi-stage to be trusted in detail. But the overall result 
- very roughly equipartition large-scale magnetic field in a few dynamical times - hinges on just 
one simple observation. Namely, we note that the number of eddy turnover times at the small-scale 
end of the Kolmogorov cascade is ^ 1. For all practical purposes, = oo, leading to an 
equipartition field even starting from a zero seed. 

To be concrete, consider a merger of two gas clouds of size ~ R. The merger shocks and stirs 
the gas for a few dynamical times ~ where V is the velocity of the merger. The resulting gas 
motion is hydrodynamic: R is much larger than the mean free path of the hydrogen atom A (this 
allows to simulate the structure formation using hydrodynamics). 

The gas thermal velocity should be of order V, because numerical experiments show significant 
heating during the first mergers (Abel et al 1998). Then the gas viscosity is ~ V\. Since the 
Reynolds number is Re = ~ y ^ 1, the gas motion is turbulent. 

The Kolmogorov cascade to small scales develops. The characteristic velocity of eddies of size 
/ is Vi ~ V{j^)^^'^. The eddy turnover time at scale / is r; ~ tdi^f'^. The Kolmo gorov cascade 
cuts off when viscosity starts to dominate, ~ r^^^, giving the cutoff scale ly ~ (-f )^''^A. The 
number of turns of the smallest eddy in the cascade is ~ ^ ~ {^Y^"^ ■ 

For example, with R ~ 10^'^ cm and A ~ lO^^cmEl, we get ~ 1000. If the small-scale magnetic 
field exponentiates after each turn of the smallest eddy, the growth factor is ~ ~ 10^'''^. This 
number is clearly too big to mean anything. All it says, is that an arbitrarily weak small scale field 
will saturate in a time <C td- 

Then the magnetic field will climb up the Kolmogorov cascade, staying in rough equipartition 
with the eddies of the covered scales. The time to climb the entire cascade is dominated by the 
last step, because td + -^j^ + + ... ~ t^. 

3. The scenario 

As we said, with the amplification factors ~ 10^'''^, it does not really matter where the seed 
field is coming from. But we think the full story develops roughly by the following scheme: 

shot noise at Debye length <5 — kinetic plasma instabilities (scales above S) — small-scale 
dynamo (scales below 1^) — dynamo (scales above /j^, climbing up the Kolmogorov cascade). 



^Here and below all length scales are physical, with the numbers loosely taken from Abel et al (1998). 
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3.1. Debye length and kinetic plasma instabilities 

The Debye length 6 and the plasma frequency ojp are defined as (5 = ( ^^^^^^^ )^^'^ ~ ^> where 
T is the gas temperature, Xe is the ionization fraction, n is the gas density. With T ~ lOOOK, 
Xe ~ 10~^, and n ~ Icm^^, we get the thermal velocity V = 3 x lO^cm/s, the Debye length 
5 3 X lO^cm, and the plasma frequency LOp ~ 10s~^. 

The mean free path of the hydrogen atom is A ~ ~ lO^^cm, where a ~ 10~^^cm^ is the 
atomic hydrogen scattering cross section. The mean free path of charged particles from Coulomb 
scattering Ac ~ ~ lO^^cm is of the same order. 

Since 5 <C A, the shocks in the merging gas develop kinetic instabilities. The growth rate of 
kinetic instabilities is ~ tOp, the length scales are ~ 6. The seed for these instabilities is provided 
by the shot noise at the Debye scale: {xen6^)~^^'^ ~ 10~^. This is a sufficient seed, because the 
number of e- foldings of kinetic instabilities can be estimated as ujpy ~ 3 x 10^. The instabilities 
are mostly electrostatic, but for our purposes, it is enough to note that magnetic fields of order ^ 
will accompany any electric field in a moving plasma. 

There is a kinetic instability (Weibel) which gives magnetic fields directly (Sagdeev 1966), but 
we don't have to worry about these details here. We simply note, that due to kinetic instabilities, 
we get a shot-noise magnetic field on scales I, very roughly of order (f )~'^''^. 



3.2. Magnetic Reynolds number and the small-scale dynamo 

The plasma conductivity is o" ~ ^m'^ iSF ~ 3 x lO^^s^^. This gives the magnetic diffusivity 
rj = ~ 3 X lO^cm^s^^, and the magnetic Reynolds number Rm = ~ 10^^. 

For Rm ^ Re, the magnetic fields exists on small length scales, below the hydrodynamic tur- 
bulence cutoff scale l^. On these length scale the hydrodynamic flow acts as a linear incompressible 
deformation. On the average, these linear deformations amplify the magnetic field 

The growth rate of the small-scale dynamo on the scales below l,y is ~ r^^^ ~ Nt^'^. The 
cutoff scale of the small-scale dynamo Im is set by the resistive damping, Nf^^ ~ V^^^ giving 
Ijn ^ Re~^^^Rm~^^'^R ~ lO^'^cm. Since the number of e-foldings of the small-scale dynamo is 
~ 10^, the shot noise seed from the kinetic dynamo will suffice to saturate the magnetic field on 
scales below li, in a time <C t^. 



^ Kulsrud, Anderson (1992) give a physical discussion, Gruzinov et al (1996) give a rigorous proof. 
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3.3. Climbing up the Kolmogorov cascade 

Assuming the small-scale dynamo provides a seed on the scale ~ li,, what should we expect on 
scales I > in a turbulence with Rm ^ Re ^1? There is no real theory for this process 0, but 
we expect that the magnetic field will climb up the Kolmogorov cascade, spreading to larger and 
larger scales. In a few dynamical times, a roughly equipartition, roughly large-scale field should 
appear. 

This behavior is dictated by the following three principles. First, magnetic energy is not 
conserved, and should spread to all available scales "by ergodicity". Second, the spread of mag- 
netic energy should be approximately local in the wavenumber space. Third, the only reasonable 
saturation level of the field amplitude is the equipartition with the flow[f|. 



4. Structure formation with magnetic field 

Assume that the first dynamo does operate. What are the consequences? The synchrotron 
emission of the first shocks is negligible. It appears that the first magnetic fields are potentially 
important for just one thing, the structure formation. In particular, they might change the initial 
mass function of the first stars and black holes. 

Obviously one cannot directly simulate the kinetic instabilities and the structure formation in 
one go, because tOp^ ~ 0.1s and ~ 3 x lO^^s. One cannot even simulate the hydrodynamic part 
of the problem, because both Re and Rm are too large. 

But we think there is a way to include the first magnetic fields in the structure formation 
simulation. One should run several simulations with Rm > Re > 1, where ">" means greater by 
as much as numerically possible. The simulations should include an initial seed field of tunable 
magnitude. If the overall scenario of §3 is correct, one will find that the simulation results are 
roughly independent of the seed field magnitude, if the seed magnitude lies in a certain interval. 

This work was supported by the David and Lucile Packard foundation. 



^Because, by the "3D Zeldovich theorem", magnetic field changes the turbulence long before it reaches the large- 
scale equipartition (Gruzinov, Diamond 1994). 

^ But we must note, that in a somewhat different context of relativistic dynamo, where we predicted "very roughly 
equipartition large-scale magnetic field in a few dynamical times", Zhang et al (2009) actually find that although the 
large-scale field indeed grows at a much faster rate than t|^^, after a few td the large-scale magnetic field stays below 
the equipartition. 
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